Tomonaga-Luttinger-liquid criticality: numerical entanglement entropy approach 
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The von Neumann entanglement entropy is studied with the density-matrix renormalization group tech- 
nique. We propose a simple approach to calculate the central charge using the entanglement entropy for one- 
dimensional (ID) quantum system. This approach is applied to a couple of quantum systems: (i) ID frustrated 
spin model and (ii) ID half-filled spinless fermions with nearest-neighbor repulsion; and, it is confirmed that the 
central charge is estimated very accurately for the both systems. Also, a new method to determine the critical 
point between TL-liquid and gapped (or ordered) phases from the proposed approach is suggested. Furthermore, 
we mention that the Tomonaga-Luttinger parameter can be obtained in a like manner as the central charge, using 
the charge-density fluctuation of a part of the ID system. 

PACS numbers: 71.10.Pm, 03.67.Mn, 75.40.Mg, 71.30.+h 



I. INTRODUCTION 

Quantum entanglement is a quantum mechanical phe- 
nomenon created by the separation of interacting quantum 
particles. This fundamental concept underpins a wide range 
of research in physics, such as quantum information, quantum 
computing, and quantum gravity, etc. When one considers 
this phenomenon in a study, it is often quantified as an entan- 
glement entropy. In the field of condensed-matter physics, the 
von Neumann entanglement entropy has been frequently used 
as an order parameter to investigate quantum phase transition 
or topological order in the quantum many-body problems. In 
particular, for one-dimensional (ID) quantum systems the en- 
tropy can be directly related to the central charge of the con- 
formal field theory (CFT). 2,3 To date, a method with this re- 
lation has been developed as a powerful tool to study the uni- 
versal properties of the Tomonaga-Luttinger (TL) liquid^, in 
combined with numerical methods such as exact diagonaliza- 
tion (ED) or density-matrix renormalization group (DMRG) 5 
techniques. 

Let us consider a quantum ID periodic system with length 
L. The von Neumann entanglement entropy of its subsystem 
with length I is given as Sl(1) = — Tr;/5; log/5;, where pi = 
TrL-ip is the reduced density matrix of the subsystem and p 
is the full density matrix of the whole system. Using the CFT, 
the entropy of the subsystem with length I for a fixed system 
length L has been derived^M 
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where c is the central charge of the associated CFT and si 
is a non-universal constant. A prime objective of using this 
formula is to estimate the central charge,-^ which provides 
definitive information concerning the universality class of 
(1 + 1) dimensional system.? This estimation has been applied 
to a variety of fermionic systems; e.g., with impurity^— , 
coupled to bosonsr^ii and under the magnetic field^ etc. Al- 
though the central charge can be also obtained by examining 
the finite-size correction of ground-state energy or the low- 
temperature behavior of specific heat,- excited states have to 
be taken into account and considerable difficulty is involved 



in the ED and DMRG calculations. But, if one uses Eq.(Q]i, it 
is only necessary to consider the ground state. 

So far, two kinds of numerical methods have been prin- 
cipally adopted to extract the value of c from Eq.(Q]i: one 
is direct fitting of Sl(1) as a function of I for a fixed sys- 
tem length and the other is scaling analysis using a relation 
AS = Sl(%) — Sl,i(^-) ~ |lnp-, considering two systems 
with lengths L and L'. However, one would encounter diffi- 
culty attributable to a finite-size effect in their practical use. 
More specifically, in the former method the equation (HJ is 
exact only in the thermodynamic limit L^oo and some am- 
biguity remains in the fitting of finite-size result^ whereas, 
in the latter method the system length to be studied must be 
strongly restricted due to the incommensurability (or frustra- 
tion) of spin and/or charge fluctuations as well as the (open- 
)shell problem in itinerant systems. Therefore, the aim of this 
paper is to propose a new method that overcomes those prob- 
lems. We first derive an efficient formula for calculating the 
central charge. And then, to ascertain the validity of the for- 
mula we apply it to a couple of ID quantum systems; (i) frus- 
trated spin model and (ii) itinerant model with strong charge 
fluctuation. We demonstrate that a very accurate estimation 
of critical point between TL-liquid and gapped (or ordered) 
states is enabled with this method. Furthermore, we suggest 
that the TL parameter can be calculated in a like manner as 
the central charge. 



II. DERIVATION OF THE FORMULA 



Let us now derive the formula for estimating the central 
charge. We first prepare Sl{j) and S(j — 1) from Eq.tQ}- 
Then, by subtracting the one from the other a simple expres- 
sion is obtained: 
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Using this expression we can easily calculate the central 
charge for a ID system with a fixed system length. The main 
advantages of Eq.© are; (i) the non-universal constant si 
does not appear explicitly, (ii) only two values of Sl (I) at the 
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TABLE I: System-size dependence of ci and C2 obtained from 
Eqs.lfJJ and <[3j for the ID S = | frustrated Heisenberg model. 
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middle of system, where the finite-size correction to Eq.© is 
the smallest^ and the DMRG calculation is the most accu- 
rate, are used, and (iii) a unique value of c is estimated for 
each system length and it enables us to perform a systematic 
extrapolation to the thermodynamic limit. Also, an alternative 
expression can be derived from the second derivative of Eq.© 
with respect to I. Since <Sx(Z) is symmetric about I = -j, i.e., 
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Note that Eq.([3j can be also derived by applying 
In [cos (£)] ~ into Eq.©. 

m. APPLICATION TO QUANTUM SYSTEMS 

A. frustrated spin chain 

In order to confirm the validity of Eqs.© and ©, we apply 
them to a couple of ID quantum systems. The first system 
considered is the S — | frustrated Heisenberg chain. The 
Hamiltonian is given by 
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where Si is a spin-i operator at site i. The parameters Ji and 
J2 are nearest-neighbor and next-nearest-neighbor antiferro- 
magnetic exchange interactions, respectively. This system 
has been extensively studied both analytically and numeri- 
cally*^— the ground state is of a dimerized zigzag-bond state 
for 0.241 < J2/J1 < 0.5 and of the Majumdar-Ghosh state 
with incommensurate spiral correlations for J2/J1 > 0.5; ac- 
cordingly, the spin gap opens when J2/J1 ^ 0.241. While 
for J2/ J\ < 0.241, the system is described as a TL liquid and 
the central charge is expected to be unity (c = 1). In recent 
years this system has been frequently used as a touchstone of 
numerical methods.^— 

By the DMRG method we study the systems with lengths 
L = 32 to 144 keeping m ~ 25L density-matrix eigenstates. 
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FIG. 1: (Color online) System-size dependence of ci as a 
function of J2/J1. The arrows indicate maximum positions 
(J2,max/Ji, ci.max). The system length is L = 32 (red), 64 (blue), 
96 (green), and 128 (orange) from top to bottom. 



Note that an appropriate ID array for the construction of the 
PBC is necessary to obtain highly-accurate result (see Fig. 1 
of Ref.l28l). In this way, the largest discarded weight is u>d ~ 
2 x 10~ 12 in the renormalization procedure and the central 
charge converges at least six digits as a function of m for a 
fixed system length. As examples, the DMRG results of ci 
and C2 at J2/J1 =0 and 0.2, where the system is in the TL- 
liquid phase, are shown in Table |I]for several system lengths. 
We can see that both of the quantities c\ and C2 converge very 
fast with increasing the system length and they can be easily 
extrapolated to 1 in the thermodynamic limit. In fact, even for 
L = 32 the largest deviations from c = 1 are only 0.5% and 
0.7% for ci and C2, respectively. A finite-size correction is 
always smaller in c\ than in C2, so that we make use only of 
ci, namely Eq.©, for the finite-size-scaling analysis below. 

Now we demonstrate how one can estimate a critical point 
between the TL- and non-TL-liquid (or gapped) phases with 
Eq.©. This is a good enduring test of Eq.© because the 
use of AS is no longer allowed due to the system-size depen- 
dent frustration of the system ©. This system belongs to the 
Gaussian universality class (c = 1) for the TL-liquid phase 
and c < 1 is expected for the gapped phase from the renor- 
malization in the massive region^ In Fig. Q] the DMRG re- 
sults of ci is plotted as a function of J2/J1 for several system 
lengths. In the renormalization group the correction of central 
charge is expressed as c = 1 + C(/ 3 ) where / is the coupling 
constant of marginal operators, 3 -'- so that a 'transition point' 
for a fixed system length may be given by a maximum of c\. 
Actually, we find that all of those curves have a maximum 
at some point (J2,max/Ji, Ci ;max ); on the left-hand side of 
the maximum point they becomes flatter and flatter, and gets 
closer to Ci = 1 with increasing the system length; whereas 
on the right-hand side, c\ drops down abruptly as a function 
of J-xj ' J\. Hence, the maximum position must approach the 
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FIG. 2: (Color online) (a) Maximum position J2,max/Ji as a func- 
tion of l/L for systems with even number of sites from L = 32 
to L — 144. The red dashed line denotes a linear least-square fit 
with a weight function w(l/L) = L. (b) Maximum height ci jmax 
as a function of l/L 2 for systems with even number of sites from 
L = 32 to L — 144. The red dashed line shows a simple linear 
fitting. 



TABLE II: System-size dependence of ci and C2 obtained from 
Eqs.l(2j and $3$ for the ID half-filled spinless model with nearest- 
neighbor repulsion. 



V/t = 

L Ci C2 



V/t = 1.9 

Cl c 2 



30 1.001119 1.002954 1.005450 1.007293 

70 1.000199 1.000535 1.002542 1.002879 

110 1.000082 1.000218 1.001818 1.001955 

150 1.000041 1.000114 1.001467 1.001540 
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FIG. 3: (Color online) System-size dependence of ci as a function of 
l/L 2 for V/t — 0. The system lengths are taken at intervals of 20 
sites from L = 50 to L = 150. The red dashed line shows a simple 
linear fitting. Inset: similar figure as the main one for V/t = 2. The 
blue dashed line is guide to the eyes only. 



critical point J2 = J2, c as the system length increases and 
reach there with c\ = 1 in the thermodynamic limit; namely, 
J2,max -> ^2,c and ci, max ->■ 1 as L -> CXJ. 

In Fig. [2ja) the maximum position J2,max/Ji 1S plotted as 
a function of l/L for systems with even number of sites from 
L = 32 to 144. Although it oscillates due to the inconsis- 
tency between the periodicity of spin wave and the system 
length, i.e., system-size dependent frustration, its amplitude 
decreases rapidly with increasing system length. Therefore, 
an extrapolation of J2,max/Ji to the thermodynamic limit is 
practicable and we obtain J% tC / J\ = 0.24112778 by linear 
least-square fit with a weight function w(l/L) = L [without 
the weighting, Ji, c / J\ = 0.2414 ± 0.001]. This value is con- 
sistent with the previous estimationsi2L22 Another thing to be 
examined is an extrapolation of the maximum height to the 
thermodynamic limit. The maximum height ci. max is plotted 
as a function of l/L 2 in Fig.|2jb). Interestingly, the data seem 
to be on a straight line at least for longer systems and we ob- 
tain c = 0.999992824 in the thermodynamic limit by a linear 
fitting. It implies that the central charge scales as 

ci = 1 + 0(1/L 2 ). (5) 

This can be interpreted as the 0(1 /L 2 ) corrections originated 
from the x = 4 irrelevant fields in the CFT.— No sooner J 2 



moves away from J 2 c than the system-size dependence of C\ 
starts to deviate from Eq.©. It means that the logarithmic cor- 
rections are completely eliminated only at this critical point. It 
is consistent with the fact that the effective Hamiltonian of the 
system is purely Gaussian only at the critical point. In this 
manner we can easily determine whether or not the logarith- 
mic corrections are present for arbitrary ID quantum system. 

B. spinless fermions with nearest-neighbor repulsion 

It would be of importance to examine whether Eqs.© and 
(O are applicable to itinerant fermion systems. Thus, as our 
second test we consider the ID half-filled spinless fermions 
with nearest-neighbor repulsion. The Hamiltonian is given by 

H = -t^2(c\cj+ h.c.) + V ^2 inrij, (6) 

where c\ (cj) is a creation (annihilation) operator of a spinless 
fermion and rij (= c|c;) is the corresponding number opera- 
tor. The repulsive interaction V (> 0) is assumed to act only 
between neighboring sites The particle density is fixed 

at n = 1/2. This system can be mapped onto the exactly solv- 
able Heisenberg XXZ chain and it is known that a transition 
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FIG. 4: (Color online) Finite-size-scaling analysis of ci as a function 
of V/t. The arrows indicate maximum positions (Vmax/t, Cl,max )• 
The system length is L — 30 (red), 70 (blue), 110 (green), and 150 
(orange) from top to bottom. 



between TL-liquid metallic and charge-density-wave insulat- 
ing phases occurs at V/t = 2(= V c /t). 

We study the systems with lengths L = 30 to 150 keeping 
m ~ 40L density-matrix eigenstates. This way, the largest 
discarded weight is ~ 8 x 10~ 12 in the renormalization 
procedure and the central charge converges at least seven dig- 
its as a function of m for a fixed system length. Tablellllshows 
the DMRG results of ci and c 2 at V/t = and 1.9, where 
the system d6} is in the TL-liquid phase, for several system 
lengths. In common with the case of the frustrated spin chain, 
the convergence of c\ and cn with the system length is very 
fast and the deviation from c = 1 is always below 1% for 
L > 30. In Fig. [3] we plot the values of c\ as a function 
of 1/L 2 for some values of V/t. At V/t — 0, we see that 
Ci is scaled by Eq.© and extrapolated to c\ = 1.00000062 
in the thermodynamic limit; however, for V/t > Eq.© is 
no longer fulfilled due to the occurrence of the logarithmic 
corrections. As an example, the finite-size scaling of c\ at 
V/t = 2, where the logarithmic corrections are maximum, is 
shown in the inset of Fig. [3] 

So, let us check the behavior of c\ on V/t near the crit- 
ical point. In Fig. [4] the values of ci is plotted as a func- 
tion of V/t for several system lengths. As seen in the frus- 
trated spin chain, each of the curves have a maximum at a 
point (Vmax, ci jm ax)- Note that the peak position is exactly 
V m ax/t = 2 for all system lengths. This is because the 'crit- 
ical point' is independent of the system length in this model 
(|6}. Quite interestingly, it means that the 'critical point' is not 
affected by the logarithmic corrections. Thus, we decide that 
Eqs.© and ® are useful for ID itinerant system as well. Sep- 
arately from this paper, it has been verified that this method is 
successfully applied to a ID spinless fermions with boson af- 
fected hopping. 34 

Lastly, we mention that the TL parameter K p , which is one 



TABLE III: System-size dependence of K P (L) obtained from 
Eqs.([8j for the ID half-filled spinless model with nearest-neighbor 
repulsion. The extrapolated values to the thermodynamic limit (L — > 
oo) and the exact values [Eq.((9jl] are also shown. 



L V = V = l.8 V = 1.9 

30 0.9926566 0.5212309 0.4940692 

50 0.9973670 0.5476091 0.5202473 

70 0.9986855 0.5587141 0.5314329 

90 0.9992341 0.5648065 0.5376011 

110 0.9995373 0.5683950 0.5417063 

L -> oo 1.0000593 0.5843504 0.5596052 

Exact 1.0000000 0.5838163 0.5562247 
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FIG. 5: (Color online) Finite-size-scaling analysis of K P (L) as a 
function of 1/L 2 for V/t = 0. The system lengths are taken at 
intervals of 20 sites from L = 30 to L = 110. The red dashed line 
shows a simple linear fitting. Inset: finite-size-scaling analysis of 
KpiV) as a function of 1/L for V/t — 1.8 and 1.9. The blue and 
green dashed lines are fitting with function K P (L) = I\ p + aL -7 . 
Through all the calculations L/2 is kept to be odd, and therefore, 
a < (see text). 



of the most notable quantity to study the TL-liquid properties, 
can be estimated in a manner similar to above. The TL param- 
eter is related to the charge-density fluctuation of subsystem 
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where A and fi are non-universal constants. The function 
Fl(1) behaves similarly to Sl(1) with /. Then, in the same 
way as Eq.© is derived, we easily obtain for L ^> 1 
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Since the 0(L 2K p) correction in Eq.(|7| oscillates on alter- 
nate sites, the set of -Fl(-j) and — 2) is a more proper 
choice than that of F L ( -§ ) and F ( § - 1 ) . The values of K p (L) 
estimated from Eq.® at V = as well as near the critical 
point are shown in Table Hill With increasing L, K p {L) seems 
to approach the exact value 



p 2arccos[-V/(2t)]' v ' 

We can easily perform the finite-size-scaling analysis with a 
fitting function K p (L) = K p + aL r (1^7<2) where a 
and 7 are fitting parameters. From Eq.®, we know a > 
(a < 0) for L/2 =even (odd) and expect a simple relation 
7/2 K p if the logarithmic corrections are small. The fit- 
ting analyses are shown in Fig. [5] As expected, 7/2 = 1 
is obtained at V/t = 0. Near the critical point, we have 
7/2 = 0.5312 (0.4982) for V/t = 1.8 (1.9) and they have 
a kind of deviation from the exact values of K p due to the 
existence of logarithmic corrections. Nonetheless, the extrap- 
olated values are reasonably close to the exact ones (see Ta- 
ble mil l. The remaining error in the thermodynamic limit is 
only about 0.1% (0.6%) for V/t = 1.8 (1.9). 

Here, we briefly comment on applying the open boundary 
conditions (OBC). Thus far the periodic boundary conditions 
are assumed; however, in general the OBC are more prefer- 
able in the DMRG study for accurate and convenient calcu- 



lations. In this regard, with applying the OBC the function 
Sl(1) oscillates like Eq.(|7]» due to the Friedel oscillation and 
the analysis with Eq.© is unreasonable. If one can suppress 
the Friedel oscillation sufficiently at the center of the system, 
for example, by using the sin 2 -deformed OBC techniquer^^ 
the central charge might be calculated by taking S(^) and 
S(j - 2) as Eq.® even with the OBC. 



IV. CONCLUSION 

A numerical approach to calculate the central charge for ID 
quantum system using the entanglement entropy is proposed. 
This approach is applied to the frustrated spin chain and the 
half-filled spinless fermions with nearest-neighbor repulsion 
to check the validity. It is confirmed that the central charge is 
estimated very accurately for both the models even with short 
system length. We also suggest a new method for determining 
the critical point between TL-liquid and gapped (or ordered) 
phases by using the proposed numerical approach. Further- 
more, we demonstrate that the TL parameter can be calculated 
in a like manner as the central charge. 

We appreciate Satoshi Ejima for enlightening discussion. 
We also thank Aroon O'Brien and Yohei Fuji for helpful dis- 
cussions. 
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